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AN IMPROVED METHOD SOP. CALCULATING THE DYNAMIC 
TESPONSE OF FLEXIBLE AIRPLANES TO GUSTS 
By Abbott A. Putnam 

SUMMARY 


A set of equations baaed, on the first mode of wing bending is 
presented for determining the dynamic responses of an airplane wing 
structure as induced by gusts. The representation of the airplane is 
such that the aerodynamic damping of the vibratory motion of the wing 
is separate from that of the vertical motion of the airplane as a 
whole. An easily evaluated solution of the equations for a unit— Jump 
forcing function is also presented which may be built up to determine 
the response to any forcing function. A chart is included by which, 
for a conventional airplane, a typical forcing function for any 
gradient gust may be chosen quichiy. 

Satisfactory agreement within the limits of accuracy of the data 
was found in a comparison of calculated response of a semirigid test 
model with experimental results. 


INTRODUCTION 


The dynamic response of airplane wings to gusts has become of 
greater interest ae the size and speed of transport and bomber type 
airplanes has increased. Eussner (reference l) sets forth the general 
equations for determining this response, but the complexity of the 
solution of these equations has led most investigators to make 
simplifying assumptions in order to reduce the equations to a form 
which is more readily solved. In reference 2, the assumptions made 
reduce the equations to those for a biplane equivalent which simulates 
the fundamental mode of bending of the wing of the actual airplane. 

A comparison of calculations with test results in reference 2 
indicates that the method used yields maximum values which are in 
good agreement with the test results but that the calculated time 
histories of the responses do not agree well. The reason for the 
disagreement was found to be that, with the biplane equivalent, the 
aerodynamic-damping coefficients of the wing bending and the airplane 
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aa a whole, which vaiy In relative Importance, are approximated by 
a single constant in each damping tern. 

In view of the limitations of the method of reference 2, the 
method presented herein was devised to provide aerodynamic damping 
that Is more representative of the actual case. Although it was 
derived independently from the concepts of the effects of stability 
and unsteady lift used in reference 2, the flexural system was 
considered in much the same way as it was considered by Williams 
and Hanson in reference 3. 

Solutions of the equations of the present method for a unit- 
jump forcing function are also presented. These solutions may be 
built up Into the response to arbitrary forcing functions with much 
less labor than when, as in reference 2, a direct solution for each 
forcing function is performed. 

Also presented is a comparison of test results for the semi- 
rigid model of reference 2 with results calculated by tbe present 
method. In addition, a chart is presented to enable a typical 
forcing function for any gust gradient distance to be chosen quickly 
for a conventional airplane. 


SYMBOLS 


t time, seoonds 

p air density, slugs per cubic foot 

dC L 

- — slope of lift curve, per radian 
da 

T forward velocity, feet per second 

S gross wing area including area Intercepted by 

fuselage, square feet 

S v net wing area, square feet 

Sf intercepted fuselage area, square feet 

y coordinate along span of wing, feet 

m(y) mass of wing per unit span, slugs per foot 

c(y) wing chord as a function of wing spanwise station, feet 
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c B root chord of wing, feet 

c.{. tip chord of wing, feet (tapered wing) 

g(y) normalised deflection curve of wing, 

6^ absolute displacement of fuselage, positive upward, feet 

5^ deflection of wing tip with respect to fuselage from steady 
flight condition, positive upward, feet 

8 absolute displacement of any station on wing, positive 

upward, feet 

M total mass of airplane, slugs 

Uf fuselage mass, slugs 

My wing mass, slugs 

My. equivalent wing mass, slugs 

X damping coefficient of airplane, pound- seconds per foot 

fuselage damping coefficient, pound-seconds per foot 

X w wing damping coefficient, pound -seconds per foot 

Xy equivalent wing damping coefficient, pound- seconds per foot 

f w natural frequency of wing, fuselage fixed, cycles per second 

f Wf natural frequency of wing-fuselage system about nodes, cycles 
per second 

K spring constant, the upward force on fuselage due to both 

wings when in assumed deflection shape, per unit tip 
deflection,, pounds per foot 

Fy fraction of forcing function acting on wing 

Ff fraction of forcing function acting on fuselage 

Ate"^ forcing function assumed for conventional airplane 

An.V 

A - at bt = 1 or A = Wbe £n_ 
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h = 

tp 

to 

1 
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arbitrary load-factor increment that the airplane would 
experience if it had no vertical motion as it traversed 
the gust, g 

total weight of airplane, poundB 


time to reach forcing function peak from entry, seconds 

time to reach gust peak from entry, seconds 

unit -.jump wing forcing function M » 0 for t < 0 
1 for t > 0 ) 


£ - 


Mkf 

XM-p 

M^v 






KM 

MfMw. 


Z • real root of equation z 3 + (t) + 0)z^ + (7 + rjp ) z + 7 *» 0 
real part of complex roots of foregoing equation 


(♦--*¥*) 


a 


' imaginary part of complex root of equation represented "by Z 


(« - /? - * 2 ) 


A = i]f - 
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r = 3Z 2 + 2( n + (3)z + (y + n£) 


p 


FfM 

M f 


FfMr) 

Q = 

E = PZ 2 + QZ + y 


F^M F f M My 

L =^~W~^ 


N 


F v Mf3 F f M 
Mw e 


B a LZ + N 


METHOD OF ANALYSIS 
Equations 

In the derivation of the present equations for the dynamic 
response of airplane wings to gusts, certain of the initial assump- 
tions are the same as those of reference 2, namely: 

(a) The loading is symmetrical 

(h) Only the first mode of bending is important 

(c) The deflection at any point on the wing may he expressed 

with the necessary accuracy as g(y) times the tip 
deflection relative to the fuselage where g(y) is 
variable with span but Independent of time 

(d) Damping varies linearly with vertical velocity 

Instead of the biplane equivalent of reference 2, however, the air- 
plane Is considered as shown in figure 1. With this representation, 
the equation for the summation of the vertical forces on the fuselage 
under the unit— Jump— type forcing function Ff J. is: 
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Mf - 


d 2 6f 


dt* 


d6f -1 

o + *-f - K 5a = Ff 1 


CD 


and the equation for the summation of the vertical forces aver the 
entire wing under the corresponding forcing function. F w J is: 


2 Ttlp 

-'root 


(2a) 


When the Indicated integrations are performed, equation (2a) reduces 
to: 




d 2 8 f 

dt 2 


+ 


dB f 

~dt 


+ 1%, 


dt 2 8 


d &d 

dt 


+ 



(2b) 


The advantage of using the present concept of the airplane as 
shown in figure 1 over the biplane equivalent of reference 2 is 
immediately apparent when the damping terms of equations (l) and 
(2b) are examined. The damping of the vertical motion of the air- 
plane as a whole is now separated from the damping of the vibratory 
motion of the wing, which thus eliminates the necessity for the 
compromise made to represent the damping when the biplane equivalent 
is used. This result Is achieved by redefinition of the mass and 
mass distribution and by basing the damping of the equivalent wing 
mass on the vibrational velocity rather than on its absolute vertical 
velooity , 


The solutions of equations (l) and (2b) under action of the 
uni tr- Jump forcing function for the fuselage acceleration increment, 
wing— tip acceleration increment relative to the fuselage increment, 
and the wing-tip displacement relative to the fuselage with the 

- d5f d5j 

initial conditions of steady flight, when * — - » — * e o are: 

dt dt 


dt 2 


) Bin Sfflt JJ- 


If; 

fl \ 


cos nt 


Pjr + PZ + Q 


+ T 


1 


( 3 ) 
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M Zt 
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M* 


[- 


Gw) «•*<* 


( 4b 'N 

Wz + ?W 


x 

sin - nt 

M 



1 


(5) 


The acceleration increment of any point on the wing relative to 
the ground, which ie of interest in certain design problems, may be 
found from the equation 

c)2g d^Sf d^s^ 

It 5 = It 5 + 8 y "at 5- (6) 


If the applied gust is assumed to be unifonn across the span of 
the airplane so that the distribution of forcing function is similar 
to that of the damping of the vertical motion, that is, if 


1 _ El - Ze 

X Xf Xw 


equations (4) and ( 5 ) may be simplified to: 


^!4-i|s 2 i zt + i n 

dt 2 M l r 


1 



^ cos 



- _ M L 

4 ' x2 r 



ijfza zn\ , x. 

•fiT + T~J sin "M nt 

+ e**Y - cos A nt 
\ M 



sin 


X ^ 

s n V 


(7) 

(8) 


If, in a given case, it is felt that torsion of the wing 
would invalidate the results, a first approximation to the response 
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of the first combined bending- torsi on mode may be obtained by the 
redefinition of the spring constant, the equivalent wing mass, and 
the wing damping coefficient to take into account the added degree 
of freedom. Terms which contain the redefined constants then merely 
replace the terms previously used when equations (3), (k) , and (5) 
are evaluated. 


Evaluation of Constants 

Damping coefficients .- The damping coefficients are evaluated 

from the relation \ => O.75 75 SV where —— 75 SV is the term 

for the steady lift per unit vertical velocity increment. An 
effective damping factor of 0.75 is inserted to take care of the lag 
in lift and variable damping of the actual case. The simplifying 
assumption of a constant effective damping factor rather than the 
variable factor of the exact theory is Justified by 4 a series of tests 
and computations reported in reference 2. The damping coefficients 
are therefore 

dCV 


VS. 


w 


0.75 


dC, 

da. 


VS 4 


and 




= 2(o.75 


dC I 

da 


0 £ 


c(y) g(y) dy 


The coefficient X 


, Wq includes the variation of wing velocity 


relative to the fuselage in accordance with the assumed wing deflec- 
tion curve. If the deflection of the wing is assumed to vary with 


square of the distance from the root, X™ will equal 

0 


elliptical wings, and 



T 


for tapered wings. 


Spring constant and equivalent wing mass .- When the equivalent 
wing mass M„. is determined from the definition 
ft ip® 

M We = 2 / m (y) g(y) dy, the spring constant K may be computed 

e '“root 

from M We and either of the natural frequencies f w or f^. 


If the 
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fuselage is fixed and no damping is present, the solution for the 

1 r k 

natural frequency of the wing leads to f «*■ = k- \ t. — . In the 


natural frequency of the wing leads to f v ^ y y— - . ' In the 
case where .the wing-fuselage combination is free to vibrate about the 
natural nodes, the solution gives the frequency ~Wf 58 if- Mvr Q Mf 


Forcing function ,— The forcing— function ratios Ff> and F w 
may be evaluated from the assumption that the forcing-function forces 

If 

are distributed proportionally to the damping forces; thus Ff = 

^w 

and F w = — and, although this assumption may not be exactly correct, 
it allows Ff and F w to remain constant throughout the time history. 

S olution for roots of auxiliary equation.— In order to 
determine the solution of the differential equations, the cubic 

z3 + (q + P)z 2 + (7 + qB)z + 7 « 0 

must be solved for the three roots, Z, + i ft, and' — ift, For 
typical cases, the value of Z haB been found to be between -1.000 
and —1.020, and can be determined with sufficient accuracy by 
straight line interpolation of the value of the cubic at 
z = —1.000 and z =-—1,010. The value of t is then — 

2 

and 0 is \ I — ~ t 2 . 


Determination of Response to Gusts 

In order to determine the response of the flexible airplane to 
a given gust, a forcing function representing the forces on the air- 
plane when considered to be rigid is combined with the response of 
the flexible airplane to the unit- Jump forcing function by use of 
Duhamel’s integral. Since equations (l) and (2b) contain terms 
representing the damping of the vertical motion of the airplane as 
a whole, the forcing function representing the gust need oniy contain 
the effects of the imposed gust and of the pitching stability of the 
airplane. Although the function may be computed for each gust shape 
and for each airplane, it was found (reference 2) that for most 
conventional airplanes the forcing function could be represented by 
a curve of the general form Ate"^. In order that, for a given 
airplane, the function represent a specific gust shape, a relation 
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must be established between it and the shape of the gust. It has 
been found that if a rigid airplane model is flown through a gust of 
a standard shape (that is, one having a linear increase of gust 
velooity to a maximum value which is maintained), reaching the 
maximum at the time tQ after entering the gust, the acceleration 

record of the airplane will show a maximum at approximately the 
time tQ (reference 4). The desired relation is, therefore, one for 

which the maximum acceleration of the flexible airplane, if considered 
rigid, occurs at tQ. The equation containing this relation 




A t f. 

o . M ^ 

2 H t F + r — 

M t F 



( 10 ) 


may be obtained from the equation for the reaction of the rigid 
airplane to the impressed forcing function 


M + X “f = Ate~ M (11) 

at 2 « 

by solving for the acceleration and determining the time tQ to 

peak acceleration. If the dimensionless valueB j~tp and *tQ are 

used, the solutions of equation (10) may be expressed as the single 
curve in figure 2, Thus, values of tp, which are the inverse of 

b, may be obtained directly from tQ. 


APPLICATION OF METHOD 


Since the absolute values of the response of an airplane to a 
known true gust velocity cannot be determined at present, the arbitrary 
multiplying constant A is used in the forcing-function equation. 

By application of the methods outlined in reference 2, however, the 
ratios of the dynamic responses to the responses under so-called 
static conditions may be computed. These ratios, which are inde- 
pendent of the arbitrary constant, are then used as multiplying 
factors to the corresponding stresses and acceleration increments 
determined by normal static design procedure to determine those 
responses under dynamic conditions. 
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COMPARISON OP THEORY AND TESTS 


As £ measure of the accuracy which might be expected from this 
theory, responses have been computed for the semirigid model reported 
in reference 2 and are compared herein with test results. The 
quality of the test results and the factors which may cause differences 
between theory and test results have been discussed in considerable 
detail in reference 2 and will not be mentioned further, 

A linear deflection curve was used for g(y), the deflection of 
the wing, and it was assumed that the lift and damping forces on the 
fuselage were negligible and could be set equal to zero. The spring 
constants for the two frequencies tested and other coefficients were 
evaluated from the original test data but may be obtained from the 
equivalent values givsn in reference 2 . 

In figure 3 are shown the theoretical curves for the variation 
with gust gradient distance of the ratio of the maximum wing tip 
deflection to the maximum fuselage acceleration for two wing 
frequencies, and the test results are plotted thereon for comparison. 

It should be noted that, the straight horizontal lines obtained are 
for the test hinged wing with no fuselage lift and damping. In the 
general case, the line would .have some curvature, the direction and 
magnitude of which would depend on all the airplane and gust-shape 
parameters , 

At the higher frequencies in figure 3 there seems to be a 
tendency for the curves to fall a little low for both wings. This 
result could indicate that the actual frequency of the wings was 
not quite so high as 26.1 oyoles per second but, because of the scatter 
of the points, nothing definite can be concluded. Consider only 
the left wing, which according to reference 2 produced the set of 
wing data more likely to be cori'ect. The points are noted to fit 
a straight horizontal line. In the case of curves for the lower 
wing frequency, the scatter appears larger but, on a percentage 
basis, is not. The tendency is for the data for the right wing to 
fall high and that for the left wing to fall low. Since it appears 
(reference 2 ) that the right wing hinge tended to buckle, the 
indicated differences would be expected. Again, no definite con- 
clusions can be drawn, although the straight horizontal line seems to 
be a good fit to the test points . 

Sample test wing— tip-deflection curves are compared with the 
calculated ourves for three gust gradient distances for a wing 
frequency of 13*5 cycles per second in figure 4, and for a wing 
frequency of 26.1 cycles per second in figure 5 * It will be noted 
that in all the curves the data are plotted as the ratio ~ 5 &/ 6 g ^ 
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so that prediction of the magnitude of the forcing function from 
the gust velocity is not necessary • 

The agreement "between, the shapes of the calculated 'and the 
experimental curves is satisfactory in every case, especially when 
it is considered that the general forcing function Ate"kt was 
used in the calculations and that many unknown factors can cause the 
test results to deviate from those for the assumed condition. The 
timing of the peaks, which depends largely on the 'natural frequency 
of the wings, is good, except for the sharp-edge gust, where the 
use of an average effective 4.6-chprd gradient distance, as suggested 
in reference 2, for computing all responses for the test model would 
affect the results . The amount of damping of the various curves, 
however, seems to "be somewhat underestimated in the computation. 

Since the fuselage-acceleration curves of the tests show even more 
damping present than the calculations, it appears that a small 
amount of fuselage damping should "be included. In the case of the 
conventional airplane, this fuselage damping has "been allowed for in 
the theory by use of the intercepted wing area, and deviations 
between theory and test from this cause should be minor. 

Comparison of the results shown in figures if- and 5 with similar 
results given in figures 13 and l4 of reference 2 indicates that the 
method of accounting for the damping distribution in the present 
paper yields calculated tiae histories of results that are more in 
conformity with actual conditions. Although the two methods yield 
about the same results for the maximum responses to single gusts, the 
time histories determined for single gusts by the present method 
are believed to be more accurate for use in building up the reactions 
to repeated gusts. __ 


CONCLUDING REMARKS 


A comparison of test results on a semirigid model with the 
results of calculations made by the method presented indicates 
that an improvement has been made in the calculation of the time 
histories of wing deflection and acceleration increments without 
impairing the accuracy of the calculation of maximum values. In 
addition, the solution of the basic equations for the unit- lump 
forcing function, together with the presentation of a chart for 
determining a suitable forcing function for a given: gust, has 
materially reduced the time necessary for making a set of 
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calculations of tho dynamic response of the vings of a conventional 
airplane to a gust. 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. January 22, 19^7 
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Figure J. ~ Comparison of calculated and experimental 
values of the ratio of maximum wing -tip def/ecfion 
ho maximum fuselage acceleration for [/arious gust 
gradient distances. 
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(a) Sharp -edge gust with 2.4- chord gradient distance. 
Figure 4. - Comparison of experimental and calcu/ated wing 
de flections due to a gust, fw P - 13.6 cps. 
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(b) Gust with 1.8- chord gradient distance. 
Figure 4.~ Continued. 
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(c) Gust with 2L2- chord gradient distance . 
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fa) Sharp -edge gust mth Ld -chord gradient di stance. 
Figure 6. - Comparison of experimental and calculated wing 
deflections due to a gust. = 13.5 cps. 
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Figure o.~ Continued. 













